Breakdown of Dynamical Scale Invariance in the Coarsening of Fractal Clusters 
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We extend a previous analysis [PRL 80, 4693 (1998)] of breakdown of dynamical scale invariance 
in the coarsening of two-dimensional DLAs (diffusion-limited aggregates) as described by the Cahn- 
Hilliard equation. Existence of a second dynamical length scale, predicted earlier, is established. 
Having measured the "solute mass" outside the cluster versus time, we obtain a third dynamical 
exponent. An auxiliary problem of the dynamics of a slender bar (that acquires a dumbbell shape) is 
considered. A simple scenario of coarsening of fractal clusters with branching structure is suggested 
that employs the dumbbell dynamics results. This scenario involves two dynamical length scales: 
the characteristic width and length of the cluster branches. The predicted dynamical exponents 
depend on the (presumably invariant) fractal dimension of the cluster skeleton. In addition, a robust 
theoretical estimate for the third dynamical exponent is obtained. Exponents found numerically are 
in reasonable agreement with these predictions. 
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I. INTRODUCTION 

Nonlinear dissipative systems, driven out of equilib- 
rium, decay to an equilibrium state after the driving 
agent is switched off or depleted. As long as a freely 
decaying nonlinear system is far from equilibrium, the 
relaxation dynamics are non-trivial and it is natural, in 
simple cases, to look for dynamical scaling and univer- 
sality. A wide class of nonlinear relaxation problems ap- 
pears in the context of phase ordering dynamics |l],|^] . In 
the present work we explore a new aspect of phase order- 
ing in systems with a conserved order parameter. This 
aspect appears when the minority phase has long-range 
correlations and represents (at least at early times) a frac- 
tal cluster (FC). Although such an initial condition does 
not result from a quench from high to low temperature (a 
standard setting of phase-ordering dynamics), it is by no 
means artificial. There are many two-phase systems that 
exhibit morphological instabilities and ramified growth 
at an early stage of their dynamics, and coarsening at a 
later stage. A canonical example is provided by diffusion 
controlled systems, such as an overcooled liquid or super- 
saturated solution. The stage of morphological instabil- 
ity and its implications in this system have been under 
scrutiny ^-^]. If some noise is present, a DLA-like FC de- 
velops at this stage 0. The subsequent surface-tension- 
driven coarsening of this FC is unavoidable if the system 
is isolated so that the total amount of available mass or 
heat is finite. Until recently, this later stage had received 
only a limited attention. Irisawa et al. |^| carried out 
Monte-Carlo simulations of diffusion-controlled coarsen- 
ing of a two-dimensional DLA cluster and found a power 



law with a non-trivial exponent for the cluster perimeter 
as a function of time. More recently the first results of 
investigation of the coarsening of two-dimensional DLA 
clusters as described by the Cahn-Hilliard (CH) equation 
were reported j^|. These results are briefly reviewed in 
the following. 

A crucial issue in the theory of phase ordering pro- 
cesses is the presence (or absence) of dynamical scale 
invariance Plf^. Dynamical scale invariance implies that 
there is a single dynamical length scale X(t) such that 
the coarsening system looks (statistically) invariant in 
time when lengths are scaled by X(t). It was found Ref. 
P| that dynamical scale invariance breaks down during 
the coarsening of DLA clusters as described by the CH- 
equation. On the other hand, the coarsening dynamics 
apparently exhibit scaling, by which we mean power laws 
in time. These power laws (with non-trivial exponents) 
were found for the cluster perimeter |^,^| , and for the dy- 
namical length scale that shows up in the Porod-law part 
of the equal-time correlation function ||. We will call 
this dynamical length scale "the first correlation length" . 
The absolute values of the dynamical exponents for the 
cluster perimeter and for the first correlation length are 
close to each other. The gyration radius of the cluster was 
found to be constant (within possible logarithmic correc- 
tions) . The mass dimension of the cluster was also found 
to be constant (on a shrinking interval of distances). The 
last two findings indicate that the mass transport is es- 
sentially local at this stage of coarsening. We conjec- 
tured in Ref. || that an additional dynamical length scale 
(with an exponent larger than that of the first correlation 
length scale) must show up in the coarsening morphol- 
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ogy. Furthermore, we speculated that the two different 
dynamical length scales are the average width and length 
of the cluster branches. 

Fractal coarsening occurs in many physical systems. 
Two-dimensional fractal fingering, observed in a Hele- 
Shaw cell with radial geometry (for a recent review see 
Rcf. pCfl ), exhibits coarsening at a late stage of the ex- 
periment, when forcing of the more viscous fluid by the 
less viscous fluid stops. Fractal coarsening has been un- 
der scrutiny in the context of sintering, in particular of 
silica aerogels [p"l|-|l3[ . Additional examples are provided 
by thermal relaxation of initially fractal grain boundaries 
|l4| and by smoothing of fractal polymer structure in the 
process of polymer collapse jl5| . It is interesting, to what 
extent fractal coarsening is universal. 

In this paper we significantly extend the analysis of 
Ref. || of the coarsening of DLA fractals as described 
by the CH equation. In Section 2 we will report nu- 
merical evidence for the existence of an additional dy- 
namical length scale and find the corresponding (second) 
dynamical exponent. We also introduce in Section 2 an 
additional measure of the coarsening dynamics: the to- 
tal "solute mass" content outside the cluster, and find 
the corresponding (third) dynamical exponent. In the 
rest of the paper we will try to develop some theoretical 
understanding of our numerical results. To this end, we 
consider, in Section 3, the coarsening dynamics of a single 
slender bar. The results of this analysis are employed in 
Section 4, where a simple scenario of coarsening of a FC 
having branching structure is suggested. In this scenario 
two different dynamical length scales are present: the 
characteristic width and length of the cluster branches. 
The corresponding dynamical exponents are calculated; 
they are found to depend on the (presumably invariant) 
fractal dimension D of the cluster "skeleton" . A robust 
theoretical estimate is also obtained for the third ("so- 
lute mass") dynamical exponent. These predictions are 
in agreement with our numerical results. A change of 
sign of the third dynamical exponent is predicted at a 
critical fractal dimension D cr n =4/3. Section 5 includes 
a summary and discussion. 

II. COARSENING OF DLA CLUSTERS: 
NUMERICAL RESULTS 

We start with a brief description of our simulations and 
diagnostics. We solved the (dimensionless) CH equation 

^ + \v 2 {V 2 u + u-u 3 ) =0 (1) 

numerically by discretizing it on the domain Q: < x < 
512 , < y < 512 with periodic boundary conditions. An 
explicit Eulcr integration scheme was used to advance 
the solution in time, and second order central differences 
to discretize the Laplace operator. With a mesh size 
Ax = Ay = 1 no preferred directions emerged in the 



computational grid, due to the truncation errors; a time 
step At = 0.05 was required for numerical stability. 

DLA clusters [|l6| (like the one shown in Fig. [I], up- 
per left), with radius of order 250, were prepared by a 
standard random-walk algorithm on a two-dimensional 
square grid and served as the initial conditions for the 
minority phase u — 1 . To prevent breakup of the clusters 
at an early stage of the coarsening process, we followed 
the technique of Irisawa et al. ^ and reinforced the ag- 
gregates by an addition of peripheral sites. Comparing 
the correlation functions before and after the reinforce- 
ment, we verified that the reinforcement did not spoil the 
fractal properties of the cluster. 

It is convenient to introduce the density field 

u(r,t) + l 
P(M) = ^ ' 

which varies between and 1. We identified the cluster 
as the locus where u(r,t) > 0, or p(r,t) > 1/2. The 
coarsening process was followed up to a time t = 5, 000. 
Typical snapshots of the coarsening process are shown 
in Fig. 0. One can see that smaller features of the FC 
are "consumed" by larger features, while the large-scale 
structure of the cluster is not affected. 




FIG. 1. Coarsening of a DLA fractal cluster as described 
by the CH-equation. The upper row corresponds to t = 
(left) and 34.7 (right), the lower row to t = 329.3 (left) and 
4,900 (right). 

To characterize the coarsening process, several quanti- 
ties were sampled and averaged over 10 initial configura- 
tions. In this paper we will analyze the following quan- 
tities: 

1. Gyration radius of the cluster. 

2. Circularly averaged equal-time pair correlation 
function, normalized at r = 0: 
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C(r,t) 



(p(r + r',t)p(r\t)) 



(2) 



3. An estimate of the cluster perimeter P(t), defined 
as the number of broken bonds between the aggre- 
gate sites. 



4. The "solute mass" outside the cluster: 



p (r, i) dx dy 



(3) 



p(r,t)<l/2 



Quantities 2 and 3 were computed after binarization of 
the phase field data: the value of p = 1 is attributed to 
all sites belonging to the cluster, while p = is attributed 
to the rest of sites. 

As we have already noticed, the gyration radius of the 
cluster remains constant (within possible logarithmic cor- 
rections) until the latest available times, so we will con- 
centrate on the rest of the measurements. 



A. Equal-time pair correlation function 

An analysis of the r-dependence of the equal-time cor- 
relation function at different times shows that coarsening 
operates only at small and intermediate distances. In the 
following we will consider separately two regions of dis- 
tances. 



the following we do not show the standard deviations of 
the fitting parameters if they are less than or equal to 
unity in the last significant digit presented. The dynam- 
ical exponent a — 0.26 differs from the Lifshitz-Slyozov 
value of 1 /3 observed in those cases when phase-ordering 
processes exhibit dynamical scale invariance |l|,p|,pj| . 
In contrast to such processes, no a priori form for the 
correction to the power-low fit for l/li(t) is available, so 
we consider the corrected-power-law fit (Q) as empirical. 
It works well though (see Fig. ||). 




FIG. 2. Dynamics of the equal-time pair correlation func- 
tion C(r, t) at small and intermediate distances for time mo- 
ments t = (a), 34.7 (b), 516.5 (c), 1992 (d) and 4900 (e). 



1. Small distances 



The linear behavior of the correlation function at small 
distances (the Porod law), clearly seen in Fig. ^ yields 
the first correlation length l\{t) and corresponding dy- 
namical exponent. In a two-phase system with a sharp 
interface, the first correlation length li is the average 
minimum distance between a randomly chosen point of 
the cluster and the interface. We will interpret l\ as the 
typical width of the cluster's branches (see Fig. [l]). 

To determine h(t), we approximated C(r,t) on the in- 
terval 0.7 < C(r,t) < 1 by a linear function 1 — r/li(t). 
Fig. H shows the inverse correlation length l/li(t) versus 
time, and a corrected-power-law fit 
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A, 



h(t) t a + B h ' 
of this dependence, with 

a = 0.26, A h = 0.44 and B h = 2.0 



(4) 



(5) 



This fit was obtained on the interval 30 < t < 2000 
that spans from the time when quasi-equilibrium sharp 
interfaces have already formed until the time when the 
system size becomes important (see below). Here and in 
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FIG. 3. The inverse first correlation length l/h(t) versus 
time (circles) and its fits (Q) (solid line) and Q (dashed line). 



A simpler alternative (chosen in Ref. B) is to use a 
pure power-law fit. For the last decade of time, 490 < 
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t < 4900, it gives 



h(t) = C h t c 



where 



a 



0.20 



and 



C h = 4.6 



(6) 



(7) 



The difference between the vaiues 0.26 and 0.20 is dis- 
appointingly large; it gives a measure of the maximum 
uncertainty of this exponent, caused by systematic er- 
rors. Similar uncertainties occur for other dynamical ex- 
ponents that we find in this work. Our conclusion that 
a is smaller than the Lifshitz-Slyozov value 1/3 is unaf- 
fected by this uncertainty. In Section 4 we will present 
a simple coarsening scenario that gives a theoretical pre- 
diction for a which agrees, for DLA clusters, with the 
value 0.26. 



2. Intermediate distances 



the last available time of our simulations, t = 4900, the 
ratio of r c /L is already about 0.3. Therefore, for more 
reliable results we should limit ourselves by t = 2000, as 
we have already done in the fitting of h(t). 




A close inspection of the correlation function C(r, t) 
at intermediate and large distances show that there are 
small (within 5 to 10%) changes there. These changes re- 
sult from small systematic variations of the cluster mass 
in the process of coarsening (see below). They make it 
difficult to perform accurate measurements at intermedi- 
ate and large distances, therefore some normalization at 
these distances is necessary. We normalized the correla- 
tion function C (r, i) at different moments of time to its 
values at r = 120: 



C(r,t) 



C{r = 120,t = 0)C(r,t) 
C{r = 120, t) 



(8) 



The results obtained by using the normalized function 
C(r, t) are not sensitive to the exact value of r chosen for 
the normalization, as long as it is large enough. 

Figure [| shows C(r,t) at different moments of time. 
The log-log plot helps to decide on the range of dis- 
tances and times we can work with. It is seen that C (r, t) 
chang function of time, only at distances smaller 

than some correlation radius r c (t) that increases with 
time. At distances r > r c (t) C(r,t) stays very close to 
its initial value C(r, 0). As it is our aim to investigate an 
intermediate asymptotic coarsening regime related to the 
fractal structure of the cluster at t = 0, we should work 
on a (shrinking) interval of scales where, at t = 0, C(r, 0) 
exhibits a power-law behavior p2|j . A reasonably accu- 
rate power-law fit can be achieved at t — on the interval 
3 < r < 150, while beyond r = L = 150 (the upper cutoff 
of the DLA cluster) finite-size effects become large. The 
fitting function is cr~ 5 , where S — 0.30 and c = 0.74. 
Correspondingly, the fractal dimension of our DLA clus- 
ters is D = 2 — <5 = 1.70, a reasonably accurate value in 
view of the relatively small size of our system. The in- 
equality r c (t) <C L puts an upper limit on the coarsening 
time that we can still work with. Fig. shows that for 



FIG. 4. Dynamics of the normalized equal-time pair cor- 
relation function C{r,t) [Eq. (^)] at time moments t = 
(a), 3.65 (b), 34.7 (c), 516.5 (d), 1992 (e) and 4900 (f). The 
dashed line shows the power-law fit 0.74 r -0 ' 30 on the interval 
3 < r < 150). 




FIG. 5. The second correlation length h(t) versus time 
(circles), and its fit (^) (solid line). 

The normalized correlation function C (r, t) exhibits an 
additional dynamical length scale that appears at inter- 
mediate distances. Recall that, at sufficiently large dis- 
tances, the tails of C{r, t) at different time moments co- 
incide with C(r, 0) (see Fig. |J). We define l 2 (t) as the 
minimum value of r\ such that for r > r± C (r, t) is less 
than C(r, 0) by no more than 10%. Such values of ^(i) 



4 



exist only for t > 800, and they are shown, on the log-log 
plot, in Fig||. As the available time interval in this case 
is quite short (less than one decade) , we are forced to use 
the last available time t = 4900 and limit ourselves to a 
pure power-law fit: 



a P = 0.20 and C P = 3.5 ■ 10 4 



(15) 



where 



l 2 (t) = C l2 1* 



(3 = 0.39 and C h = 1.5 



(9) 



(10) 



Though the standard deviation is small, the short time 
interval does not guarantee a high precision of the ex- 
ponent 0.39. A big difference between the exponents a 
and P is, however, beyond doubt. We will interpret ^(t) 
as the characteristic branch length which serves as the 
time-dependent lower cutoff of the fractal "skeleton" of 
the coarsening cluster. 

Equations (0) and (|l(]) show an interesting relationship 
between l\ and li\ 



hit cx t° 



(11) 



The exponent 0.98 is close to unity, and we will return 
to this observation in Section 4. 

It should be noted that normalization of C(r, t) at large 
r, which helped us to extract the second dynamical length 
scale, is no more artificial than the widely used normal- 
ization at r = 0. The only real need for any normaliza- 
tion of the equal-time pair correlation function is non- 
constancy of the cluster mass in time. When investigat- 
ing the dynamics at small distances, it is convenient to 
normalize the correlation function at r = 0. When inves- 
tigating the dynamics at large distances, it is convenient 
to normalize C(r, t) somewhere in the tail. 



B. Cluster Perimeter 

Figure || shows the dynamics of the cluster perimeter 
P(t). A decrease of the perimeter of a cluster under con- 
dition of (approximate) conservation of the cluster mass 
is a clear manifestation of coarsening. The corrected- 
power-law fit of P(t) is 



P(t) 



A, 



t a ? + Bp ' 



(12) 



where 



a P = 0.26, A P = 7.0 -W 4 and B P = 2.0. (13) 
A pure power-law fit for the last decade of time gives 



P(t)=C P r ap , 



(14) 



where 



The absolute values of the exponents ap and a coincide 
for the same type of fit. The inverse slope of the Porod- 
law-part of the correlation function should indeed scale 
as the perimeter, because each element of interface which 
size is much larger than the distance r contributes inde- 
pendently to the correlation function M . 
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FIG. 6. Cluster perimeter P(t) versus time (circles) and 
its fits (Q (solid line) and (0) (dashed line). 



C. "Solute" Mass Outside the Cluster 

Let us use, for a moment, the language of physics of 
liquid solutions. As all of the solute mass at t = is 
concentrated in the DLA cluster, rapid dissolution of the 
solute from the cluster edge occurs first. The cluster 
mass will start to decrease with time. Unless the area 
fraction of the "fractal phase" is too small the dissolu- 
tion stops, and the dissolved material precipitates back 
on the coarsening cluster [Q . The mass of the dissolved 
material decreases at this stage, while the cluster mass 
slowly increases, asymptotically approaching a constant 
value. The late stage of this regime should be qualita- 
tively similar to Ostwald ripening, where the dynamics 
of the dissolved material is responsible for a correction to 
the Lifshitz-Slyozov scaling behavior jlTj . 

These arguments give a qualitative explanation to our 
numerical results on the solute mass outside the cluster 
versus time, M s (t). Fig. |^ shows a log- log plot of M s (t) at 
sufficiently late times (the rapid "dissolution" observed 
at earlier times is not shown). The same figure shows 
two fitting functions to M s (t) . A corrected-power-law fit 
is 



M s {t) = 



A 



Ms 



^ + B M . 



(16) 



where 
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7 = 0.24 , A Me = 3.64 • 10 4 and B Ms = 2.4 . (17) 
A pure power-law fit for 500 < t < 4900 is 

M s {t)=C Ma t- 1 , (18) 

where 

7 = 0.16 and C Ms = (1-45 ± 0.02) • 10 4 . (19) 
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dynamical scale invariance is caused by the effective lo- 
cality of mass transfer which manifests itself in the (ap- 
proximate) conservation of both the mass, and gyration 
radius of the cluster in the process of coarsening. 

Looking for a scenario with broken dynamical scale in- 
variance one should, first of all, identify the nature of the 
second dynamical length scale. Snapshots of the coars- 
ening process (Fig. [|) indicate that it might be the av- 
erage length of the cluster branches Having made 
this assumption, we should verify that the typical width 
and length of the cluster branches indeed show different 
dynamical scalings, the length growing faster than the 
width. 

An important element of the scenario that we want 
to explore is the shrinking dynamics of a single slender 
bar of phase u = 1 ("solid"), evolving under the CH- 
equation ([l]) in the "liquid" phase u = — 1. Special sim- 
ulations show (see below) that, in the process of shrink- 
ing along its main axis, the bar acquires the shape of a 
dumbbell, and the "balls" at the ends of the dumbbell 
expand with time. In this Section we will derive scal- 
ing relations for the time-dependent parameters of the 
shrinking dumbbell-shaped bar, and then compare them 
with numerical simulations. 



A. Slender bar dynamics: theoretical estimates 



FIG. 7. Mass of "solute" outside the cluster M s versus 
time (circles) and its fits ( |l6| ) (solid line) and ( |l8| ) (dashed 
line) . 

One can expect that the fractal coarsening dynamics 
look simplest (that is, corrections to dynamical scaling 
are small) when the cluster mass is already almost con- 
stant. This condition is only weakly satisfied in our sim- 
ulations. For example, even at t = 2000 the solute mass 
outside the cluster reaches about 10% of the cluster mass. 
This may be the reason for the relatively large uncertain- 
ties of the numerical values of the dynamical exponents 
that we have found. Another possible reason is the fi- 
nite size effects related to a moderate fractal range of the 
DLA-clusters that we used in simulations. 



III. SLENDER BAR DYNAMICS 

In the rest of the paper we will try to get some quan- 
titative understanding of our numerical results. As no 
satisfactory theory for coarsening of FCs is available, 
one can try to formulate a simplified scenario and com- 
pare its predictions with simulations. The simplest pos- 
sible scenario assumes dynamical scale invariance, that is 
the presence of a single relevant dynamical length scale 
[^3|,0. We have seen that, in the case of CH-dynamics, 
this scenario disagrees with simulations. Breakdown of 



Our approach to the slender bar dynamics employs 
a modified version of the rigorous asymptotic sharp- 
interface theory developed for the CH-equation P,p5|. 
This theory (which can be called "Laplacian coarsening" ) 
requires all characteristic length scales in the problem to 
be much larger than the domain wall width (which is of 
order unity), but much smaller than the characteristic 
diffusion length Id ~ t 1 / 2 . Under these assumptions, 

• the local normal velocity of the moving interface 
of the bar is equal to the difference between the 
normal components of the gradient of the density 
field p(r, t) outside and inside the interface, 

• the density field is represented by two harmonic 
functions: one outside, the other inside the bar, 

• there is a Gibbs-Thomson matching condition at 
the moving interface and no-flux condition at the 
external boundary. 

This formulation enforces exact conservation of the bar 
area in the process of coarsening. On the other hand, the 
condition of an infinite diffusion length is too restrictive. 
It is not satisfied in any of the state-of-the-art numerical 
simulations of phase ordering as described by the Cahn- 
Hilliard equation. Fortunately, it is sufficient to require 
in practice that the diffusion length be larger than the 
characteristic coarsening length(s) of the problem, and 
this is what we shall do. 
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We assume that the bar length is much larger than 
its width and use the aspect ratio of the bar as a large 
parameter. Limiting ourselves to order-of-magnitude es- 
timates, we will assume that the density field inside the 
bar has already settled down, so that we can treat the 
bar interior simply as a region where p(r,t) = 1. 

What is the physical picture of the bar dynamics? Be- 
cause of the large aspect ratio of the bar, the density gra- 
dient will be largest near the dumbbell ends and small 
elsewhere [see Eq. (|20| ) below]. The solute, released from 
the bar ends, will therefore be transported to the liquid, 
and the bar will be shrinking along its main axis. Let us 
assume that the total mass of the solute outside the bar is 
relatively small (see below). Then most of the dissolved 
solute will be reabsorbed not far from the bar ends. This 
effectively local dynamics is the reason that the bar ac- 
quires the shape of a dumbbell. The "balls" at the ends 
of the shrinking dumbbell will effectively travel along the 
main axis of the bar like snowballs, accumulating mate- 
rial along their motion and growing in size. 

Going over to a quantitative analysis, we consider a 
half-infinite bar and denote the width of its planar part 
by A. We assume that, despite the fact that the ball at 
the end of the dumbbell may have a complicated shape, 
it can still be characterized by a single time-dependent 
length scale R(t). (Recall that what we are after is 
order-of-magnitude estimates, rather than a complete so- 
lution.) Let the bar be placed along the x-axis as shown 
in Fig. |8|. In this Figure xo(t) is the time-dependent po- 
sition of the bar edge (that is, the length reduction of the 
bar in the process of its shrinking). Introduce the polar 
coordinates r and <j> with the origin at the bar end. 




FIG. 8. Setting for the dumbbell dynamics. Dashed lines 
are the lines of flow calculated from Eq. (po[). 

The density field p must be a harmonic function that 
vanishes at the planar part of the dumbbell interface 
(that is, at (j> — > and <f> — ► 2ir) and is of order \/R(t) 



at the "ball" interface (for definiteness, at cf> = ±ir/2). 
Far enough from the bar end, r 3> R(t), the bar can be 
considered as a thin sheet [we assume that A <C R(t)] ■ 
The solution of the Laplace equation for a thin sheet is 
elementary [B7J , and we can estimate the density as 



p(r, <f>,t) ~ C[R(t)r}- 1/2 sim>/2), 



(20) 



where C is a constant of order unity. Employing the 
boundary condition at the ball interface, we have ex- 
tended the thin-sheet approximation to the limit of its 
applicability, but this can only affect the value of the 
constant C. Dashed lines in Fig. || show the streamlines 
of the mass flow obtained by taking the gradient of the 
density field @. 

Calculating the total mass of the "solute" with Eq. 
(f20|), one can see that it diverges at large r. This diver- 
gence arises because one replaces a diffusion equation by 
the simpler Laplace equation in the asymptotic theory of 
the CH-equation §J|5|. The diver gence can be cured by 
introducing an upper cutoff R m ax in Eq. (p0[). The cut- 
off is the smallest of the two lengths: the diffusion length 
Id ~ t 1 / 2 and the system size Lq. 

There are two important consequences of Eq. J2C|). 
First, the density gradient |V/?| is indeed largest near 
the bar end. As the result, the ball emits material from 
its edge and retreats. Second, the total mass flux to the 
planar part of the bar is finite, as it is proportional to 
r~ 3 / 2 and therefore converges at large r along the bar 
edge. The second property implies that most of the so- 
lute, emitted by the bar end, is reabsorbed by the less 
curved part of the "ball" . Again, this requires that the 
solute mass outside of the bar is small, and we will return 
to this condition later. 

Now we proceed to obtain simple scaling relations for 
the parameters of the shrinking half-infinite bar. Em- 
ploying the (approximate) constancy of the bar mass, we 
have 

x (t)A ~ [R(t)} 2 . (21) 

The mass flux out of the ball can be estimated as 

m~ \Vp\ r „ R R~ [Rit)}- 1 . (22) 

As a result of the outflow and reabsorption of this mate- 
rial the ball travels along the x-axis. The characteristic 
time tr it takes the ball to pass the distance comparable 
to its size is tr ~ m/m ~ [R(t)) 3 , so the ball speed is 



i ( t) „ M „ mr 2 

Equations @ and (§|) yield 

x (t)x (t)A ~ 1 , 

which follows 

xo(t) ~ A~ 1/2 1 1 ! 2 . 



(23) 

(24) 
(25) 
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Then Eq. Q2J) yields 



(26) 



Therefore, shrinking of a slender bar in the CH- 
equation exhibits dynamical scaling with exponents 1/2 
(for the bar length) and 1/4 (for the ball size). In the 
following we will verify Eq. ( p5[ ) numerically and then 
employ it in our fractal coarsening scenario. 

Let us return to the condition of the relative smallncss 
of the solute mass outside the bar. Using Eq. (|2(]) with 
an upper cutoff, we can estimate the solute mass outside 
the bar as [R{t)]- 1 / 2 R%l x . This should be much less 
that the ball mass: 



[R(t)]-'/ 2 RV 2 ax ^[R(t)¥ 



(27) 



At not too large times, t < Lq, it is the diffusion length 
Id ~ t 1 / 2 that should be taken for the upper cutoff R ma x- 
In this regime the condition (E71) can be rewritten as 



A > t 1/5 . 

At later times, t > Lq, we should put i?„ 
the condition (^?j) becomes 

A»^ 2/5 t^. 



(28) 
Lq, and 



(29) 



The criteria fl28|) and (|29|) coincide (by order of magni- 
tude) in the intermediate regime t ~ Lq and yield 



A >L 



2/5 



(30) 



Eq. ( |30| ) is the most stringent criterion, so we should re- 
quire that it holds. On the other hand, we should assume 
that an initially rectangular bar has already developed 
its dumbbell shape. This requires R(t) ^> A. Using Eq. 
(prf), we can rewrite it as A <C t 1 ^ . Overall, the dumb- 
bell scalings (E5J) and (E6[) are expected to hold when 



L 2 / 5 « A « *V3 



(31) 



B. Slender bar dynamics: numerical simulation 

The dynamics of a slender bar were simulated with the 
same CH equation ([!]). We started with a rectangular- 
shaped bar (u = 1 inside the bar, u = — 1 outside). The 
bar sizes were 512x8, and the bar was placed in the center 
of a 1024 x 256 rectangular box. The CH-equation (|l|) 
was solved with periodic boundary conditions. Only one 
quadrant was actually simulated because of symmetry 
with respect to the x and y axes. 

Snapshots of the bar evolution at time moments t = 0, 
8154, 26810 and 88180 are shown in Fig. §. They con- 
firm the "dumbbell picture" qualitatively. One can see 
that most of the bar keeps its planar shape and constant 



width. Balls are formed at the ends of the bar and travel 
along the bar's main axis growing in size. 
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FIG. 9. Dumbbell formation and dynamics at t = 0, 8154, 
26810 and 88180 (from top to bottom). 

Fig. |l0| shows the numerically obtained time- 
dependence of the bar length reduction Xo(t). A pure 
power-law asymptotics of this quantity, predicted by Eq. 
(p5|), is expected to show up at very late times, when the 
typical size of the balls R(t) becomes much larger than 
A. At intermediate times we can expect a correction of 
order A/R(t). For A = 8 used in our simulations, this 
correction is significant until the latest available times 
~ 10 5 . Taking it into account, we can fit the numerical 
data for xo(t) by the function 



x Q (t) = C\t 1 / 2 + C 2 t 1 ^ 



(32) 



This fit (with C\ — 0.6 and C2 = —1.6) and its leading 
term are shown separately in Fig. [Io[ and a good agree- 
ment is observed. Reintroducing the A-dependences, 
predicted by the theory, we can rewrite Eq. (p2|) as 



x (t) = 1.7 (i/A) 1/2 - 0.95 (A • t) 1 ' 4 + 



(33) 



In the following Section we will use only the leading term 
ofEq. (|3|). 

The parameters chosen for this simulation made it pos- 
sible to reach the dumbbell scaling regime. Indeed, the 
right inequality in Eq. ([H]) is satisfied for any of the 
time moments shown in Fig. |^, except t = 0. Taking the 
system length Lq — 124 (the distance between the planar 
part of the bar and the boundary), we see that the left 
inequality in Eq. ( |3l| ) requires A»7. The results of 
the bar dynamics simulation for A = 8 indicate that a 
usual (not strong) inequality is sufficient. We observed 
breakup of the bar into fragments for A = 4, when this 
condition is violated. 
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FIG. 10. The bar length reduction xo(t) versus time (solid 
line); fit ( |32| ) (dashed line) and its leading term (dotted line). 

The results of our investigation of the slender bar dy- 
namics are interesting in their own right. One should 
work with domains of a simple shape in order to un- 
derstand the basics of coarsening dynamics. It is a sin- 
gle spherical droplet of the minority phase that usually 
serves as a "test object" for coarsening dynamics 
Our results show that a slender bar represents an instruc- 
tive alternative. First, the slender bar dynamics exhibit 
locality of mass transport. Second, two different dynam- 
ical length scales appear. Neither of these two properties 
is present in the single droplet dynamics. 

In conclusion of this section we cite earlier papers pq ] 
where analytical solutions for the Laplacian growth prob- 
lem were obtained in a slender bar geometry. Differences 
in motivation (diffusion-limited growth versus diffusion- 
limited coarsening) lead to important differences between 
these models and our model. First, the boundary con- 
dition far from the bar in Refs. |2g] corresponded to a 
non-zero flux of material, while this flux is zero in our 
model. Second, papers |2(| did not account for surface 
tension at the bar interface. 



IV. "FRACTAL SKELETON" SCENARIO OF 
COARSENING 

As it is evident from Fig. ([!]), a two-dimensional DLA 
cluster preserves, in the process of coarsening, its branch- 
ing structure. We will characterize the typical cluster 
branch at time t by a width a(t) and length b(t). Later on 
we will identify a(t) and b(t) with the two time-dependent 
correlation distances li(t) and hit), respectively. We 
are interested in the coarsening regime when the total 
mass (area) of the cluster is almost constant. Define 
the skeleton of the cluster (at fixed time) by tending all 
branches widths to zero. Coarsening of a FC in our sce- 



nario involves disappearance of the shortest branches of 
the cluster and an increase of the width of the remaining 
branches so that the cluster mass remains constant. We 
assume that the cluster skeleton, rather than the cluster 
itself, preserves its fractal structure, with the same frac- 
tal dimension as at t = 0, in the process of coarsening 
p8[ . The lower cutoff of the fractal skeleton is the typical 
branch length, and it grows with time. The upper cut- 
off remains constant. Another important assumption is 
that each individual daughter branch evolves like a single 
slender bar of the previous Section, until the time when 
its aspect ratio becomes of order unity. Then this branch 
rapidly shrinks and disappears, "injecting" its material 
into the parental branch. 

Now Eq. ( ^5|) can be interpreted in the following way. 
Because of its shortening, the branch with the width cto 
and length bo has a lifetime of order 

t ~ a b 2 o . (34) 

Obviously, only those branches which lifetime is larger 
than t will survive by time t. Therefore, at time t the 
typical remaining branches satisfy the scaling relation 

a(t)b 2 {t)~t. (35) 

The same scaling relation apparently holds for the two 
correlation lengths l\{t) and h{t) [see Eq. (|ll|)] . This 
indicates that a(t) and b(t) can indeed be identified with 
li(t) and h(t), respectively. 

The total mass of the cluster can be estimated as 
ab (L/b) D , and this quantity must be constant and equal 
to the initial mass M. Combining it with Eq. (|35|), we 
arrive at the following scaling relations: 

2 

o(t)~ii(t)~ {jnj t— , (36) 

Kt)~h(t)~(—) t~- (37) 
For the cluster perimeter we obtain 

f L\ D P-l 2D D-l 

p(t)~bl-\ ~M~—L—t-—. (38) 

The different dynamical exponents obtained for ait) and 
b(t) explain breakdown of dynamical scale invariance. 
Notice that the absolute values of the exponents for P(t) 
and ait) coincide. In the limiting case of D = 2 all the ex- 
ponents coincide to give the Lifshitz-Slyozov value 1/3, so 
that dynamical scale invariance is restored. This limiting 
case describes a quench through the critical point, where 
convoluted percolating interfaces (for which D = 2) are 
observed 

In addition, we can give a robust description of the 
scaling behavior of the "solute" mass outside the cluster. 
The "injection" processes outlined above causes undula- 
tions of the interface of the parent branches (clearly seen 
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in Fig. [l]). In their turn, diffusion of the "solute" and 
surface tension tend to erase the undulations. It is natu- 
ral to assume that, by time t, undulations with wave- 
lengths smaller than ~ i 1 / 3 have already disappeared 
p9| . Hence, at time t the typical wavelength (or cur- 
vature radius) of the branch undulation is of order t 1 ^ 3 . 
[Notice that, for D < 2, this curvature radius grows with 
time slower than b(t).] The "solute" density p around 
the curved interface is of order i -1 / 3 due to the Lapla- 
cian screening of "bays" by "capes". Solute with this 
density will be found within a distance of the order of 
the diffusion length Id oc t 1 ^ 2 from the cluster interface. 
Outside of this region the density will very small. The 
area of this region can be estimated as 

S d ~ l\ (J^) oc t^ . (39) 

Multiplying this quantity by the "solute" density t^ 1 / 3 , 
we obtain a dynamical scaling relation for the "solute" 
mass outside the cluster: 

4-3D 

M s oc t~ — . (40) 

Notice that for D = 2 the Lifshitz-Slyozov scaling i -1 / 3 
for the solute mass fl7j| is recovered. A striking predic- 
tion of Eq. ( poj ) is the change of sign of the dynamical 
exponent at D = 4/3. For D < 4/3 the "solute" mass 
should continue increasing with time until the time when 
the dissolving fractal degrades and the "fractal skeleton" 
model becomes inapplicable. It should be noted that es- 
timate ([It)]) is quite robust, as it is independent of most 
of the assumptions of the "fractal skeleton" scenario. For 
example, it does not use the values of the first two dy- 
namical exponents and only assumes that they are less 
than 1/2 (that is, the length and width of the branches 
grow in time slower than the diffusion length Id). 

Now we are in a position to compare the predictions 
of the "fractal skeleton" scenario with our simulations of 
DLA coarsening. This comparison is made in Table 1 for 
D = 1.70. The agreement is quite reasonable, in view of 
the uncertainty range of the exponents found numerically 
and simplicity of the scenario. 



TABLE I. Dynamical exponents: scenario vs. simulation 



Quantity 


exponent 
from scenario 


exponent 
from corrected 
power-law fit 


exponent 
from pure 
power-law fit 


h 


0.26 


0.26 


0.20 


h 


0.37 




0.39 


P 


-0.26 


-0.26 


-0.20 


M s 


-0.18 


-0.24 


-0.16 



V. SUMMARY AND DISCUSSION 

We investigated the bulk-diffusion-controlled coarsen- 
ing of DLA fractal clusters as described by the CH- 
equation. We observed that long-ranged correlations, in- 
troduced by the FC at t = 0, define a new intermediate 
asymptotics (in the terminology of Ref . J3(]] ) in the coars- 
ening dynamics: "fractal coarsening" . In order to reach 
this asymptotics, one should wait long enough so that 
quasi-equilibrium domain walls are already formed and 
mass transfer between the cluster and the majority phase 
is already weak (the cluster mass is approximately con- 
stant). On the other hand, this stage is limited at large 
times by finite-size effects [a finite value of l2(t)/L]. This 
intermediate asymptotic stage is quite long (our coarsen- 
ing scenario predicts that its duration is of order L D+1 ), 
and it certainly deserves attention. 

The main result of this work (and of the preceding Let- 
ter ||) is breakdown of dynamical scale invariance, that 
is the presence of several dynamical length scales, during 
the bulk-diffusion-limited fractal coarsening. We identi- 
fied two dynamical length scales from the evolution of the 
equal-time pair correlation function, and a third dynam- 
ical length scale from the evolution of the "solute" mass 
outside the cluster. Breakdown of dynamical scale in- 
variance is caused by the effective locality of mass trans- 
fer which manifests itself in the simultaneous (approxi- 
mate) conservation of the cluster mass and gyration ra- 
dius in the process of coarsening. Scale-invariant mass- 
preserving coarsening would obviously require shrinking 
of the FC |llj , while no shrinking is observed in this sys- 
tem. 

Looking for a simple scenario of coarsening with a bro- 
ken dynamical scale invariance, we investigated an auxil- 
iary problem of the dynamics of a single dumbbell-shaped 
domain and found "unusual" dynamical exponents 1/2 
and 1/4. Locality of mass transfer is present already in 
the single-dumbbell dynamics. We suggested a simple 
scenario of fractal coarsening in diffusion-controlled sys- 
tems with a conserved order parameter. It postulates a 
fractal skeleton with an invariable fractal dimension and 
employs the dumbbell model for the dynamics of individ- 
ual branches. In addition, a robust estimate of the "so- 
lute mass" exponent is obtained, and a qualitative change 
in the "solute mass" dynamics is predicted at D = 4/3. 
Theoretical predictions are in reasonable agreement with 
numerical simulations. 

Much more work is needed, however, before a 
more complete understanding of the fractal coarsening 
emerges. A moderate fractal range of the DLA realiza- 
tions that we worked with made it difficult to obtain 
sharp estimates of the dynamical exponents. More ex- 
tensive numerical simulations would increase the scaling 
range and test the values of the dynamical exponents. 

How can one put the "fractal skeleton" scenario of 
coarsening under additional tests? The scenario gives 
very definite predictions of the dynamical exponents in 
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terms of the (invariable) fractal dimension of the clus- 
ter skeleton. In general, fractal dimension is only one 
characteristics of a geometrical set. Therefore, one direct 
test would involve simulations of coarsening of a different 
random FC with branching structure that has the same 
fractal dimension as DLA. Further tests would involve 
simulations with random FCs having branching struc- 
tures with different (tunable) fractal dimensions. Simula- 
tions of this type were performed in Ref. |l3| for the case 
of edge-diffusion-controlled fractal coarsening. In our 
case of bulk-diffusion-controllcd coarsening these tests 
will check, in particular, the prediction of a qualitative 
change in the "solute mass" dynamics at D cr n = 4/3. 

Finally, a comparative investigation of different mech- 
anisms of fractal coarsening (edge diffusion versus bulk 
diffusion, local conservation versus global conservation, 
etc.) is needed if we want to address the question about 
possible universality classes of fractal coarsening. 
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